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ABSTRACT 

We examine the clustering of galaxies around a sample of 20 luminous 
low redshift (z £ 0.30) quasars observed with the Wide Field Camera-2 on 
the Hubble Space Telescope. The HST resolution makes possible galaxy 
identification brighter than V = 23.5 and as close as 2" to the quasar. We 
find a significant enhancement of galaxies within a projected separation of 
£ 100 fo _1 kpc of the quasars. If we model the qso/galaxy correlation function 
as a power law with a slope given by the galaxy/galaxy correlation function, 
we find that the ratio of the qso/galaxy to galaxy/galaxy correlation functions 
is 3.8 ± 0.8. The galaxy counts within r < 15 ft, _1 kpc of the quasars are too 
high for the density profile to have an appreciable core radius ( £ 100 /i~ 1 kpc). 
Our results reinforce the idea that low redshift quasars are located preferentially 
in groups of 10-20 galaxies rather than in rich clusters. We see no significant 
difference in the clustering amplitudes derived from radio-loud and radio-quiet 
subsamples. 

Subject headings: quasars; galaxy clustering 
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1. INTRODUCTION 



Over the last two decades, it has been well established that quasars are associated 
with enhancements in the galaxy distribution. Historically, this provided the first direct 
observational evidence that quasars were indeed cosmological in origin (Bahcall, Schmidt, 
& Gunn 1969; Bahcall & Bahcall 1970; Gunn 1971; Stockton 1978) . Over the years, 
considerable evidence has accumulated that low redshift, z £ 0.4, quasars reside in small 
to moderate groups of galaxies rather than in rich clusters (cf. Hartwick & Schade 1990, 
Bahcall & Chokshi 1991, and references therein). At higher redshifts, there is a marked 
difference in the environments of radio-loud and radio-quiet quasars (Yee & Green 1987). 
At redshifts z £ 0.6, radio-loud quasars are often found in association with rich clusters 
(Abell richness R > 1) while radio-quiet quasars appear to remain in smaller groups, or 
perhaps in the outer regions of clusters (Boyle et al. 1988; Yee 1990). 

The galaxy environment around quasars provides many important clues as to what 
triggers and sustains their central engines. As first suggested by Toomre & Toomre (1972), 
mergers and interactions of galaxies can provide an efficient mechanism for transporting 
gas into the inner regions of a galaxy or quasar. There have been attempts to model the 
interaction/merger rates of ordinary galaxies in order to explain the luminosity function 
of quasars (De Robertis 1985; Roos 1985; Carlberg 1990) and the rapid evolution of the 
merger/interaction rate in clusters with redshift may provide a natural explanation of the 
strong evolution of clustering observed around radio-loud quasars (Stocke & Perrenod 1981; 
Roos 1981). Knowledge of the galaxy environment around quasars is also important for 
understanding the large scale distribution of quasars and how it relates to the structure 
seen in galaxy surveys (Bahcall & Chokshi 1991). 

In this Letter, we examine the galaxy environment around twenty nearby (z £ 0.3) 
bright quasars observed with the Wide Field and Planetary Camera-2 (WFPC2) of the 
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Hubble Space Telescope (HST). These fields were imaged as part of an ongoing investigation 
into the nature of the host environment of quasars (Bahcall, Kirhakos, & Schneider 1994, 
1995a, 1995b, 1996a). The exceptional resolution of the HST images allows companion 
galaxies to be detected at very close projected separations, in some cases r ~ 3 /i _1 kpc 
( £ 2"), and galaxy/star separation to be performed down to V ~ 23.5. The goal of the 
work presented here is to quantify the excess of galaxies associated with these quasars. The 
outline of the Letter is as follows. A brief description of the quasar sample used in our 
analysis is given in § 2. In § 3.1 we argue that galaxy counts are inconsistent with being 
drawn from a uniform background. In § 3.2, we strengthen this conclusion by correcting 
the counts for the background contamination. We also present the excess galaxy counts 
above the background in annuli of projected separation. From these counts, we quantify the 
amplitude of the galaxy clustering around the quasars in § 3.3 in terms of a quasar-galaxy 
spatial cross correlation amplitude. We discuss our results and their relation to previous 
work in § 4. 



2. DATA 



The sample of objects analyzed in this Letter consists of 20 of the intrinsically most 
luminous (M v < -22.9, for H = lOOkms^Mpc -1 , fi = 1) nearby (z < 0.30) radio-quiet 
and radio-loud quasars selected from the Veron-Cetty & Veron (1993) catalog. Table 1 
lists the individual quasars and their redshifts. The quasars span the redshift range 
0.086 < z < 0.29 with a median redshift of z me d = 0.18. Details of the observations and the 
resulting images for eleven of the twenty fields have been presented in Bahcall, Kirhakos, 
& Schneider (1994, 1995a, 1995b); the remaining observations will be presented in Bahcall 
et al. (1996b, in preparation). Six of the quasars (denoted with an asterisk in Table 1) are 
radio- loud, while the remainder are radio-quiet (Kellermann 1989). 
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Each quasar field was imaged with the WFPC2 through the F606W filter, which is 
similar, but slightly redder than the V bandpass (A = 5940A, A A = 1500 A). The quasars 
were positioned within 4" ± 1.2" of the center of Wide Field Camera CCD 3 (WF3). 
Simultaneous images were obtained in the adjacent CCD chips 2 and 4 (WF2 and WF4 
respectively) which together with WF3 form a "L" shaped image (see figure 1). Each 
chip has 800 x 800 pixels and an image scale of 0.0996" pixel -1 at the chip's center; this 
corresponds to spatial resolution of 2.1 ft, _1 kpc/arcsec at a redshift z = 0.20 (Qo = 1). The 
effective areas (areas not covered by pyramid shadows) of WF2, WF3, and WF4 are 1.59, 
1.60, and 1.59 sq arcmin respectively. More detailed information on the WFPC2 and its 
photometric system can be found in Burrows (1994), Trauger et al. (1994), and Holtzman 
et al. (1995a,b). The relatively long exposures (1100 or 1400 seconds), combined with the 
excellent spatial resolution, allowed galaxies to be identified in the images down to limiting 
magnitude m(F606W) < 24.5 and close as £j 2" from the central quasar. We performed 
aperture photometry on the field galaxies; circular apertures with radii of 0.3" to 10" were 
used as appropriate. 



3. GALAXY COUNTS AROUND QUASARS 

3.1. Raw Counts: Evidence for a Strong Enhancement 

If galaxies are distributed around low redshift quasars with a power law distribution, 
£(r) ~ (r/10 /i -1 Mpc) -1 ' 77 (as suggested by, e.g., Yee & Green 1987), then there will be a 
very strong enhancement of the counts within 30" (r ~ 60 ft, _1 kpc at z = 0.2) of the quasar. 
Moreover, because the centers of WF2 and WF4 are offset from the quasar, there will be an 
enhancement of galaxies in WF3 relative to WF2 and WF4. The background galaxy counts 
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rise steeply with magnitude and this will dilute any signal of excess galaxies. Much of the 
background contamination can be removed simply by counting only those galaxies with 
apparent magnitudes in the range that is likely to be physically associated with the quasars. 
A good compromise between eliminating too many associated galaxies and minimizing the 
effect of foreground/background interlopers, is to count galaxies in each field, i, which are 
in the magnitude range m*(zj) — 1 to m*(zi) + 2, where m*(zj) is the apparent magnitude 
of an galaxy at the redshift of the quasar, Z{. The mean (median) value of m*(zi) 
for our sample is 18.2 (18.1) ; the total range of m*(zi) is 16.4 < m*(zj) < 19.2. In 
computing m*(z), we have taken an galaxy to correspond to an absolute magnitude 
M*(F606W^) = —20.75 and have used the K-corrections between absolute and apparent 
magnitudes given in Fukugita, Shimasaku, & Ichikawa (1995). 

Figure 1 shows the positions of all galaxies in our twenty fields with 
Tn*{zi) — 1 < m < m*(zi) + 2; geometric corrections were applied according to 
Holtzman et al. (1995a). In the panel containing the quasar (WF3, lower left), there is 
a clear excess of galaxies (50 galaxies) relative to WF2 (upper left) and WF4 (18 and 17 
galaxies respectively). The hypothesis that the counts in all three chips are drawn from 
an underlying Poisson distribution with a common mean leads to a maximum likelihood 
mean per chip of 28.4 and a chi-squared of x 2 — 24.5; the probability that x 2 f° r one 
degree of freedom is this large by chance is only P = 7 x 10~ 7 . Thus, without any detailed 
modeling, we can rule out the possibility that the counts are random fluctuations in the 
background distribution. In the following sections, we attempt to make this conclusion 
progressively more quantitative by first modeling the background galaxy distribution and 
then introducing a model for the spatial distribution of galaxies around the quasars. 



3.2. Correcting for the Background Galaxies 



-7- 



In order to further quantify the excess of galaxy counts around the quasars, we need 
an estimate of the contribution from background galaxies. In figure 2, we show the galaxy 
counts ( d 2 N/dmdf2 , per square arcsec) versus magnitude for the "off quasar" chips WF2 
and WF4, in bins Am(FQ06W) = 0.5. The counts at m(FQOQW) > 21.5 mag are well 
approximated by a power law, log 10 d 2 N/dmdf2 = —10.8 + 0.33 m(FQOQW); the counts at 
brighter magnitudes are in excess of those obtained by extrapolating the faint power law 
fit. Because the separation of the centers of chips 2 and 4 is not large (101" corresponds to 
213 /i _1 kpc at z = 0.2), the counts in these chips are a combination of both background 
galaxies and the (relatively bright) galaxies physically associated with the quasars, and 
hence yield an overestimate of the true background. We have compared the power-law fit 
with the counts derived from the HST Medium Deep Survey (MDS) (Griffiths et al. 1994; 
S. Casertano, private communication) which covers a much larger area of the sky. The 
agreement is good, and in the remainder of this Letter, we will adopt the power law in 
figure 2 as our estimate of d 2 N/dmdfi . The agreement with the MDS is also a useful 
consistency check for systematic errors in our derived magnitudes. 

In the twenty fields, there are 11 galaxies in the range m*(zi) — 1 < m < m*(zi) + 2 
within a projected separation less than 25 ft, _1 kpc of the quasar; the total number expected 
from the d 2 N/dmdfi power law relation is only 0.99. The probability of the observed 
counts being a Poisson realization of the background is extremely small, P = 8 x 10~ 9 , 
two order of magnitudes smaller than our previous estimate obtained by neglecting the 
background. This is also a much stronger constraint than the upper limit given in Bahcall, 
Kirhakos, & Schneider (1995a), P = 2 x 10~ 2 , based on eight fields and an (over) estimate 
of the background obtained from the counts in WF2 and WF4. 

We have also computed counts within projected separations of 5 and 10 fa _1 kpc. In 
the twenty fields, we find 2/0.039 and 5/. 16 (number/number expected in the background) 
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with r < 5 ft, _1 kpc and r < 10 ft, _1 kpc respectively; the corresponding random probabilities 
are P = 7 x 1CT 4 and P = 7 x lCT 7 . There is a suggestive difference in the counts 
derived from the radio-loud and radio-quiet subsamples. The radio-quiet quasars (14 fields) 
had 0/.027 (P = 0.97) and 2/. 11 (P = 5 x 10~ 3 ) for the 5 and 10 h^kpc counts; the 
corresponding numbers for the radio-loud sample (6 fields) were 2/. 012 (P = 8 x 10~ 5 ) and 
3/.05 (P = 2 x 10~ 5 ). 

In figure 3, we show the fractional excess of galaxies above the background, 
(SN/N(r)) = (Ni(r) / Nb ti (r)) — 1, obtained by averaging the counts in the twenty fields in 
bins of 15 /i _1 kpc projected separation. Here iVj(r) is the actual count of galaxies with 
projected separations between r — Ar/2 and r + Ar/2 in the i th field and N bji (r) is the 
expected background contribution. Again, we have only considered those galaxies in the 
apparent magnitude range m^Zi) — 1 to m^Zi) + 2. 

From figure 3, we can immediately draw two conclusions. First, there is a significant 
excess of galaxies within projected separations of r < /i~ x kpc from the quasars. Second, 
there appears to be no significant difference in the galaxy counts for the radio-loud and 
radio-quiet subsamples for 10 /i^kpc < r < 100 /i^kpc. In the next section, we quantify 
the observed clustering in terms of the spatial quasar/galaxy cross-correlation function. 

3.3. Estimating the Spatial Clustering Amplitude 

A detailed derivation of the relation between angular counts and a spatial distribution 
of galaxies in terms of a cross-correlation function is given in Longair & Seldner (1979). 
Briefly, one assumes that the galaxy distribution (above the background) around the 
quasar is specified by a quasar/galaxy cross-correlation function, (, qg (r, z). The observed 
excess in projected separation in the i th field is then obtained by integrating £, qg (r, z) 
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over the redshift distribution of galaxies (dj\f/dz) i in the apparent magnitude range 
[m*(zi) — 1 < m < m*(zi) + 2], 

N-M/N, (r) - 1 - Jd6dZ (^ij^j (1) 
l[ ]l hA ] Jd9dz(dAf/dz) t ■ [i) 

In this equation, s is the comoving separation between between the line of sight, x(z), and 
the quasar located at coordinate x(z i ) ) s 2 = {(x(zi) — x(z)) 2 /G[x(z)]} + x(z) 2 9 2 ; 9 is the 
angle between the line of sight and the quasar and G(x) is a function which describes the 
degree of spatial curvature G(x) = 1 — (H x/c) 2 (£l — 1) (cf. Peebles 1980, eqn. 56.13). 
The range of integration in 9 corresponds to the angles which span the bin of projected 
separation, i.e., (r — Ar/2)/d a (zi) to (r + Ar/2)/d a (zi) where d a (zi) is the angular diameter 
distance to the quasar. 

The redshift distribution {dH/dz) i (assuming no galaxies are created or destroyed) is 
given by 

d t 

(dN/dz\ = fa) - dV , (2) 

where (f>(zi) is the integral of the luminosity function over the absolute magnitudes 
[M m (zi) — 1 < M < M*(zi) + 2] corresponding the the apparent magnitude range 
[m*(zi) — 1 < m < m*(zi) + 2], 

= / dM$(M) . (3) 

JM,{zi)-l 

The cosmological model (£1 ,A,H ) enters in the above equations implicitly in the 
volume element, angular distance, and luminosity distance. We adopt a canonical model 
with Vl = 1, H — 100 kms _1 Mpc _1 , and A = 0. The basic results of our study are largely 
independent of this choice. The redshift distribution is computed using a non-evolving 
Schechter (1976) luminosity function with a faint end slope of a = 0.97 (Loveday et 
al. 1992). The K-corrections necessary for computing the relation between absolute and 
apparent magnitudes have been taken from Fukugita, Shimasaku, & Ichikawa (1995). In 
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the results that follow, the H dependence is explicitly indicated. The derived clustering 
amplitudes remain within the stated errors as Qq is varied from 0.1 to 1.0. 

The cross-correlation function, £ qg (r,z), is usually assumed to evolve with redshift. 
A convenient model for this evolution is to assume that on small scales the clustering is 
constant in physical coordinates, i.e. the number of excess galaxies around the quasar, 
n(z)£ qg (r,z), is a constant (n(z) is average number density at the epoch z). Thus for any 
assumed shape of the correlation function, F(r), £, qg (r,z) evolves as 

1 



C q9 (r,z) 



l + z 



(4) 



The factor of (1 + z) 3 compensates for the change in the mean number density, while 
the redshift factor in the argument of F(r) is simply a matter of convention; one usually 
specifies the shape in terms of physical separations while in flU), we have specified £ qg (r, z) 
in terms of comoving separation. 

We consider three different models for F(r). The first model is a power law, 
F(r) = Bggr -1 ' 77 , with a slope equal to that of the galaxy/galaxy correlation function 
(Davis & Peebles 1983) and the amplitude taken as a free parameter. The second model is 
an exponential surface density of galaxies, fi(r) = fi exp(— r/r c ) which, after deprojection 
by a standard Abel inversion, corresponds to F(r) = fi /irr c K (r /r c ) (K (x) being the 
modified Bessel function). This model was proposed by Merrifield & Kent (1989) as a 
typical galaxy profile around centrally dominant cluster galaxies. We adopt their best 
estimate of a core radius of r c = 100 /i _1 kpc, and vary the amplitude hq. Last, we consider 
a modified Hubble profile F(r) = A [1 + (r/r c ) 2 ]~ 3 / 2 (Binney & Tremaine 1987, eqn. 2-37) 
with r c = 100 ft, -1 kpc and A a free parameter. These three models have different behavior 
at small r. Both the power law and exponential models diverge as r — >• (although the later 
does so only weakly, K (x) ~ — ln(x)), while the modified Hubble model asymptotically 
approaches a constant. 
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We varied the amplitude of the models to achieve a maximum likelihood fit to the 
excess galaxy counts in bins of projected separation of 15 /i _1 kpc. The limited number of 
fields prevented us from varying more than one parameter per model. Figure 3 shows the 
resulting best-fit models; the amplitudes of the power-law model are given in Table 2. At 
separations r >, 30 /i _1 kpc, the excess counts are relatively flat and all three models for the 
correlation function fit the data. However, the counts in the innermost bin r < 15 /i _1 kpc 
lie significantly above the counts at larger separations; the rise in the counts at small radii 
is particularly striking in the five radio loud fields. The power-law is the only model for 
£,qg( r ) considered which rises steeply enough to account for this excess. 



4. DISCUSSION 

If the quasars were distributed like typical galaxies, then the derived value of (B qg ) 
would be equal to the amplitude of the galaxy/galaxy correlation function, (B gg ) ~ 20 
(Davis & Peebles 1983). A higher value of (B qg ) suggests that the quasars lie preferentially in 
regions of above average galaxy density. Following Bahcall & Chokshi (1991), we can convert 
our values of B qg into an estimate of the typical richness of quasar galaxy environment. 
Here, we define richness as the number of galaxies associated with the quasar; this 
is given by a simple integral over the correlation function, N = 47m* J t; qg (r)r 2 dr, where 
the limits of integration are r = to 1.5 /i -1 Mpc (the traditional Abell radius), and 

10 2 h 3 Mpc 3 is the number density of galaxies. Using the values of (B qg ) in 
Table 2, we find the quasars typically reside in groups of 16-25 galaxies. These numbers 
should be compared to the typical richness of Abell clusters which have 30-49 and 50-79 
members for richness classes R = and R = 1 respectively. Moreover, this estimate is most 
likely an upper limit since there is evidence that the galaxy profile around quasars falls off 
more steeply than r~ L77 for r £ 0.25 /i -1 Mpc (Ellingson, Yee, & Green 1991). In order to 
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test the robustness of our results to a steepening of the galaxy profile, we fitted a double 
power law model with slope of —1.77 for r < 0.25 h^Mpc and —3 for r > 0.25 /i -1 Mpc. 
The best fit amplitudes, B qg , for this model increased by about 10% (still well within the 
quoted 1 — a errors), yet, the number of infered bright galaxies associated with the quasars 
decreased to 8 due to the steeper profile at large separations. 

The derived amplitudes, (B qg ), in the pure power law case are somewhat larger than 
the estimates by Yee & Green (1987). They examined the clustering of 9 radio-loud and 16 
radio-quiet quasars in the redshift range 0.15 < z < 0.30. On scales of ~ 20-500 /i _1 kpc, 
they estimated (B qg ) ps 60 ± 20 and (B qg ) ^ 42 ± 14 for radio-loud and radio-quiet quasars 
respectively. Hayman (1990) derived galaxy counts around low redshift (z < 0.3) quasars 
from the Palomar Sky Survey prints. He found that the ratio of the quasar/galaxy and 
galaxy/galaxy angular correlation functions was 3.1 ± 0.6; if the quasars and galaxies 
have similar selection functions, this translates to an estimate similar to Yee & Green 
of {B qg ) ~ 61 ± 12. French & Gunn (1983) analysed a sample 25 low-redshift quasars 
(z ^ 0.35) selected from 1.2-m Palomar Schmidt plates and concluded (B qg ) = 25 ± 12; 
they also analysed the data set of Stockton (1978, 27 quasars with redshifts z < 0.45 
selected from the red Sky Survey prints) and derived, via the same analysis, B qg = 79 ± 40. 
Our measurements are, with the exception of Yee & Green's value for radio-quiet quasars, 
consistent within the quoted errors. The slightly higher clustering amplitude we derive 
for the radio-quiet subsample may be a result of our sample being the subset of the most 
luminous quasars. 

Yee & Green (1987) found that the clustering amplitude of galaxies around radio-loud 
quasars increased by a factor ~ 3 between z ~ 0.4-0.6 and at z ~ 0.6 radio-loud quasars 
are found in environments as rich as Abell class R — 1. Optical quasars do not evolve 
as rapidly (Boyle et al. 1988), perhaps indicating a different formation scenario. It has 
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been suggested the quasars and active galactic nuclei may be triggered by interactions 
(e.g., Toomre & Toomre 1972; Stocke & Perrenod 1981; Roos 1981, Yee 1987). This offers 
a simple explanation for why the quasars are typically not found in rich clusters at low 
redshifts; the high velocity dispersion of such clusters leads to a low interaction rate. 

The HST WFPC2 is an excellent instrument for extending the present analysis to 
fainter, low-redshift, quasars. This extension would improve the counting statistics while 
also providing information regarding possible correlations of the quasar environment with 
luminosity. The imaging of of moderate redshift quasars (z £ 0.6) could be accomplished 
by HST with single orbit exposures. This imaging would provide a more direct comparison 
with previous ground based work and would increase our knowledge of the evolutionary 
history of the quasar environment. An increased knowledge of the quasar environment 
would be useful in using quasar observations to probe large scale structure at higher 
redshifts. 
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TABLE 1 

QUASAR SAMPLE 



Name 




Date 




Exp. (sec) 


My 




z 


PG 0052+251 


05 


DEC 


94 


1400 


-23 








155 


PHL 909 


17 


OCT 


94 


1400 


-22 


9 





171 


NAB 0205+02 


26 


OCT 


94 


1400 


-23 








155 


0316-346* 


20 


NOV 


94 


1400 


-24 


5 





265 


PG 0923+201 


23 


MAR 


95 


1400 


-23 








190 


PG 0953+414 


03 


FEB 


94 


1100 


-24 


1 





239 


PKS 1004+13* 


26 


FEB 


95 


1400 


-24 


3 





240 


PG 1012+00 


25 


FEB 


95 


1400 


-23 








185 


HE 1029-1401 


06 


FEB 


95 


1400 


-23 


2 





086 


PG 1116+215 


08 


FEB 


94 


1100 


-23 


7 





177 


PG 1202+281 


08 


FEB 


94 


1100 


-23 








165 


3C 273* 


05 


JUN 


94 


1100 


-25 


7 





158 
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TABLE 2 

Clustering Amplitudes 



Sample {B qg ^ 1 - a Range (B qg )/{B gg ) 

All quasars 75 60-93 3.8 ± 0.8 

Radio Loud 84 57-117 4.2 ±1.5 

Radio Quiet 72 53-92 3.6 ±1.0 
tin units of [ /i^Mpc] 1 - 77 
tt(S flP ) = 19.8( /i^Mpc) 1 - 77 
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FIGURE CAPTIONS 

Fig. I. — Mosaic of the galaxy counts in WF2, WF3, and WF4 (upper left, lower left, 
and right panels respectively) for all twenty fields. In a given field i, we count only those 
galaxies in the magnitude range m*(z;) - 1 to m^i) + 2 where m*{zi) is the apparent 
magnitude of an galaxy at the redshift of the quasar, Z; L . 

Fig. 2. — Number counts, d 2 N/dmdf2 , versus m(F606W). The filled circles represent 
the counts determined directly from the off chips in each QSO field. Poisson errors are 
shown. The line is a least squares fit to the faint counts with m(FQOQW) > 21.5. 

Fig. 3. — Average galaxy counts (5N/N(r)) of the full quasar sample and radio and 
radio-quiet subsets. The solid points are the counts computed using a power law fit to the 
faint counts of chips 2 and 4 for the background (cf. figure 2); the open symbols (shifted 
.025 in the log) show the counts derived using the background counts from the Medium 
Deep Survey. The error bars show the scatter between the different fields. The curves are 
the predicted correlation functions if the galaxy/quasar cross correlation function is a power 
law with index 7 = 1.77 (solid), an exponential surface density (dashed), or a modified 
Hubble profile (dotted). All calculations assume f2 = 1 and A = 0. 



